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ABSTRACT 
1 3 

The J Q NMR chemical shifts of graphite intercalation 
compounds have been calculated. For acceptor types, the shifts 
come mainly from the paramagnetic (Ramsey) intra-atomic terms. 
They are related to the gross features of the two-dimensional 
band structures. The calculated anisotropy is about -140 ppm 
and is independent of the finer details such as charge transfer. 
For donor types, the carbon 2pTr orbitals are spin-polarized 
because of mixing with metal conduction electrons, thus there 
is an additional dipolar contribution which may be correlated 
with the electronic specific heat. The general agreement with 
experimental data is satisfactory. 



A wide variety of molecular and metallic species may be 

intercalated between the carbon layers of graphite. Because 

of the high electrical conductivities and the nearly two- 

dimensional characteristics, these graphite i .tercalation 

1 2 

compounds (GIC) have been extensively studied. ' From the 

!3 C nuclar magnetic resonance (NMR) spectra, the chemical shifts (CS) 

may be used as microscopic probes for the electronic structures. 

The method is also universally applicable to GIC of both the 

acceptor (with molecular species such as Brg, HNO^, AsF^, 

SbCl^, AlCl^, etc.) and donor (with metallic species such as K, 

Cs, etc.) types. The pristine graphite itself has a layered 

structure with relatively weak interplanar interactions. Following 

the standard notations of the coordinate systems , the planes of the 

carbon atoms are parallel to the xy-planes and the z-axis is a 

threefold symmetry axis. The 1 ^C CS-tensor S is expected to 

be both diagonal and axially symmetric (-S is often denoted as 

4-7 

the chemical shift shielding parameter o r in the literature) . ' 

The diagonal components S__, and % will be abbreviated 

XX jj 

as S„, L. and S . Since £ = £ , there are only two independent 
z a y x y 

components, 5„ and 8 , which may be measured experimentally with 
the magnetic field B parallel and perpendicular to the z-axis. 

They will be expressed in ppm (10**^) with respect to tetramethylsilane 
(TMS) as the reference. It is convenient to define S. =( £ +2C )/3 

1 Z X 

and S = as the isotropic and anisotropic parts of the CS. 

While structural changes do occur upon intercalation, the carbon 
layer structures remain intact for GIC, hence we also expect their CS 
tensor to have approximate axial symmetry. 
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There are experimental difficulties of NMR due to the 

low natural abundances (l£) of *^C isotopes, unfavorable relaxation 

times, etc. Despite these difficulties, the general features 

8-l6 

of the experimental CS are now known. " Because large single 

crystals of graphite are not readily available, many GIC samples 

are prepared by intercalation into highly oriented pyrolytic 

graphite (HOPG) which is polycrystalline with the z-axis of 

individual crystallites in good mutual alignment throughout the 

specimen ("2° spread), though the x- and y-axes, in the graphite 

planes, are disordered from crystallite to crystallite. For 

these GIC samples derived from HOPG, sharp resonance lines are 

usually observed when B is perpendicular to the z-axis. Although 

. . 2 

structural changes do occur during intercalation , these experi- 
mental evidences (S = S„) indicate that the axial symmetry of 

x y 

the CS tensor about the z-axis is indeed retained. Also, nearly 
universal values of S ~ 40 and S vl80 have been observed for 

Z X 

all GIC of acceptor type (GICA)j these shifts are nearly independent 
of the GICA composition. On the contrary, both S_ and £ vary 
with the compositions of the GIC of donor types (GICD) . 

In the present work, we have calculated the CS for GIC 
and have considered the "band" and the "dipolar" contributions. 

We will show that the former term depends only on the gross 
features of the electronic 2ptf* and 2pir band structures and 
is not very sensitive toward the finer details such as charge 
transfers or Fermi levels. We will then correlate the "dipolar" 
term with the electronic specific heat of GIC. It will then 
be shown that only "band" contribution may be present for GICA 
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(since their electronic specific heats are very small) but both 

terms may be important for GICD. 

We will first consider the "band" contribution. The main 

contribution to CS comes from the paramagnetic (Ramsey) 

intra-atomic term while the diamagnetic term (Lamb) may be 

neglected because it is both small and nearly isotropic. For 

molecules i the Ramsey term of CS is due to the deviations 

from spherical symmetry of the electronic charge distribution ’ 

and may be calculated from the interaction between the occupied 

and empty molecular orbitals. For solids, the molecular 

. 17-19 

orbitals are replaced by Bloch functions, thus we have 

S z =(2eV/m 2 c 2 )<r- 3 > 2p Z V [Z c <VV‘V+ v N+ c ><*clU+v>J 

( 1 ) 

where the factor of 2 is included for up and down spins, c is 

light velocity, ^ is Planck's constant divided by 2tt , 

e and m are electronic charge and mass, L is the z-component 

£ 

of angular momentum operator in units of k , ^ r "^2p the 
expectation value of the inverse cube distance for the carbon 
2p-orbitals, ^ and E are the wave functions and their energies 


The subscripts v and c refer to valence (occupied) and conduction 

(unoccupied or excited) states. The functions ^ y and »^ c inside 

the square bracket have the same wave vector k. refers 

to summation over the discrete excited states. £ y refers 

to integration or average over the k-space or Brillouin zone 

(BZ). Following Strong etal. 1 ^, a reasonable choice of <r‘ ' 3> 2P 

is 1.8a = 12.2 for carbon 2p-orbitals where a^ is the 

o o 

2 2 2 2 / -3 

Bohr radius of the hydrogen atom, hence (2e * /m c )\r ^ 2 p = 

5200 ppm when E c *E y is expressed in units of eV . The expression 

for £ is rather similar with L_ replaced by L . 

Usually, informations are rather incomplete about the 

excited states, thus the energy denominators in (1) are often 

replaced by an average excitation energy dE so that the 

4 5 

closure relation may be used. %J However, the choice of 4E 
is somewhat arbitrary. For anisotrcpi.- solids such as 
graphite, there would be even more uncertainty in choosing 
two arbitrary values for and £ . 

Since CS would depend only on the general features of the 
energy bands, it may be suffice to use the tight-binding 
approximation for two-dimensional graphite where the wave 
functions of the valence state and the first excited state 
have been given by Lomer. This limited information about 
the valence state vf y , E y and the first excited state 
(subscript c=l ) » E^ (‘the lowest unoccupied state that 

satisfies the condition <+vl L = 0) permits us to give 
an over-estimate and under-estimate of CSi 
(a) We may replace the enejrgy denominator E -E in (1) by 

C V 

the common value ^E^=E 1 -E y . Using the closure relation, 
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we geti 

S z =5200Z v UE 1 )-Vf v |L 2 2 |t v > (2) 

We will refer to (2) as the average energy (AE) approximation. 
Clearly, AE is an over-estimate since is the lowest 
possible value of E c -E v and all terms in the square 
bracket of (1) are positive. The equation is similar 
(with L z replaced by L x ) for S x » However, a different 
value of AE^ must be used. 

(b) We may truncate the series in (1) by including ^ only 
under i 

S z =5200Z v <4E i r 1 |<+ v |L z |f 1 >| 2 (3) 

We will refer to (3) as the first excited state (FE) 
approximation. Clearly, FE is an under-estimate since 
all terms in (1) are positive. 

Because of the reflection symmetry in the xy-plane of 
graphite, the 7T-states are based entirely on the atomic 2p 

z 

(denoted as p Q , where the subscript 0 refers to the azimuthal 
quantum number) functions of the carbon atom, whereas the 
<T-states are based entirely on the atomic 2s, 2p + ^ and 2p^ 
(abbreviated as s, p + , p_) functions. Except for a small 
region near the center F of the BZ, the states are arranged 
from lower to higher energies in the order of O', ~n (unprimed 
denotes valence or occupied orbitals), 7 r' , a' (primed denotes 
conduction or unoccupied orbitals). The selection rules 

are t 7 t' > 7r<r' for & x and tc' for S z , thus the denominator 
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is smaller in the former case. From either (2) or (3). 
it follows that $ >S or £ <0 which is in accordance with 

a u a 

the experimental results for GIC of acceptor type. 

At the center f of the BZ, there is no mixing between the 
s and p-f unctions.*^ The pure s-f unctions are omitted here 
because they do not contribute to CS. The occupied states 
are o* ± = (f» At + Pg x )/2^ (twofold degnerate) and 
It =(P A0 + P B0^ 2 ^ where subscript A and B refer to the two 
atoms in the unit cell, t and 0 refer to the aaimuthal quantum 
number m= ±1, 0 of the atomic 2p-f unctions. Similarly, the 
unoccupied states are <r ± % = (p A -p fi± )/2^ (twofold degenerate) 
and 7r* = (p AQ -p B0 )/2^, The atoms A and B hava the same CS , 
hence it is suffice to use L zA and L^. For L^, the only 
nonzero matrix elements are the diagonal elements 

^ p A±I L zaI p A±^ ' thus the °"+ 0%' and excitations 

would contribute to S 2 . In the AE approximation, we have 

S z =(2)(5200)<T t |L 2A 2 | J - + >/(E,'-E T )=5300/(E <r t -E T ) (I*) 

On the other hand, the FE approximation is smaller by a factor of 2, 

S 2 =(2)(5200)|<’<r + | L zA | t + *>| 2 /(E_ r '-E,)=2600/(E. r *-E r ) (5) 

Similarly, the and 7TT ±' excitations would contribute 

to S'.* In the AE approximation, we nave 

S x =(5200)[2(E„'-E <r )‘ 1 <T4L xA 2 |r + > + -E„.) _1 <'^ | n- >] 

=2600 [(E pr '-£ ff .)” 2 +(E y , -E Tr )” 1 J (6) 

Again, the FE approximation is smaller by a factor of 2, 

5200)[2(EV-E T )- 1 |<r.|L xA | T -?| 2 -2(E T '.E lr )- 1 |<HL xA | x/>| 2 J 

=1300 [(E/-E,.)- 1 *(E T '-E,r 1 J (7) 
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Energy band calculations * give E ' - E- = 12 eV, 

<r w 

E ir" E <r = V"V 15 eV. For center f. we get ^=433, f x *347, 

J = =86 in the AE approximation and C =216, S =174, 

CL « A ™ A 

S a =42 in the FE approximation. 

We will now calculate CS at one of the comers P of the 
BZ. We have the same tt and 7 r' wave functions as before, 
and E w = E^' . There are now two <r -functions. One of them, 
denoted as <r^ , is twofold degenerate and is a mixture of s 
and p-f unctions. The other, denoted as o^, is a pure 
p-function. Similarly, there are the unoccupied ' and 
g* 2 ' states. In the order of increasing energies, the O’ and cr' 
functions and their energies are^i 

E l=*( E s + Ep)-if <E s- E p) 2+ 3 6H sp 2 ^- as A -fp B _. 

E 2 =E p + 3 H pp' * 

Ei • =4(E s H« p ) *iC( E s -E p ) 2 -36H sp 2 ji . «1-^B-«P A+ . ?S 4 -«P B . 

V = V 3 V’ ff 2' = ( PA--PB + )/ 2i - 

where E g and are one-electron Hamiltonian integrals for the 

2s and 2p atomic functions, and H ' are Hamiltonian 

sp pp 

integrals between neighboring atoms. They are denoted as 
E(s), E(p), H(sp) and H ( pp ) by Lomer.^ These parameters are 
all negative. The mixing coefficients * and »3 between the 2s 
and 2p functions are given by« 

- -3V 6, H ( W* [<V E p )2t 3 6H sp 2 ]*}/ 2C 
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( 8 ) 


where the normalization constant C is chosen such that u 2 +P?=l . 

The second functions of <T^ and <r^' do not contribute to 
CS of atom A since there is no angular momentum associated with 
the atomic s A functions. The and exc *' ta ' t i oris wou ld 

contribute to S z « The , ttOj ' and Tr<r 2 ' excitations 

would contribute to S x » In the AE approximation, we get 

J Z =(5200) [(E t --Ej ) - |L Z 2 |(T 1 > + (E 2 - -e 2 ) - 1 <<r 2 \L z 2 \a- 2 >] 

= (2600)[2p 2 (E 1 , -E 1 )‘ 1 +(E 2 , -E 2 )" 1 J (10) 

? x =(5200) [(e^-e, -(v E 2 rl <'^l^ 2 k 2 > t ( E 1 -V' 1 <"l L x 2 i ,r > 

= (2600)[p 2 (E ]t -E 1 )' 1 +4(E tr -E 2 )‘ 1 +(E 1 , -E ir )‘ i J ( 11 ) 

In the FE approximation, we get 

J z =(5200) [(Ej •-E 1 )''|<r<T, |L Z '>| 2 + (E 2 '-E 2 ) ■ 1 k= r 2 | L z»°' 2 ’>| 2 ] 

=5200oC 2 p 2 (E 1 ' -E 1 ) " 1 +1300 (E 2 ' -E 2 ) -1 (12) 

S x =(5200)[(E;-E 1 )- 1 |<^|L x |Tr'>| 2 

♦(Ei'-E x )- 1 |<a|L x |T 1 , >| 2 *(E 2 '-E it )- 1 \<n | L x |<T 2 *>| 2 

= (1300)[p 2 (E ji -E 1 )' 1 +i(E 7r -E 2 )' 1 + ol 2 (E 1 , -E p )' 1 +i(E 2 , -E Ir )‘ 1 . (13) 

For corner P, energy Dand calculations 20 ’ 21 gives E^- 22 , 

E 2 =-20, E w =E w> '=-9, Ej'= 5 and E 2 '=l 6 eV. In addition, it is 
necessary to evaluate the mixing coefficients <* and ^ . An 
estimate of Ep =-6 may be obtained by averaging the c r and & ' 
energies at the center T of BZ. From the values of Ej. and E^ ' , 
we get E g =-ll and H g p=-4.4. Substitution into ( 8 ) gives 
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o(=0.?6 and ^=0.65 which are consistent with the estimate 
<x 2 =|3 2 =0.5 or <* = ^*0.71 of Zunger 21 . From (10) to (13)* we get 
S =189, $ =426, $=-2. 37 for the AE approximation and £ =82, 

S’ =177 » S =-95 for the FE approximation. 

X CL 

For the symmetry point Q at the midpoint of the edge 
of BZ, a similar calculation gives £=242, £ =446, £ =-204 for the 

u X cl 

AE approximation and £,=94 , S =163, S =-69 for the FE approxi- 

mation. On averaging over the entire BZ, we get the average 

values £ =-170 and -70 for the AE and FE approximations respectively. 

3 . 

For GICA, the net effect of the charge transfer is the partial 

removal (usually several per cent) of the valence ir-electrons 

near th° corner P of the BZ. This charge transfer has no effect 

on £ since contributions to £ comes from the <rcr' excitations 
z z 

which are independent of the TT-bands. The effect on S x is also 
small. With the partial removal of the 7T-electrons , there is 
less 7 rcr' excitation. However, the vr-band is now partially 
vacant, hence some c rir excitations are also possible. There is 
considerable cancellation between these two contributions since 
the tt and tt' bands are close together near the corner P. Similarly, 
the charge transfer effects on both £, and 9 may also be small 
for GICD. 

In summary, the "band" contributions to £ and £* come 

from the entire BZ and are independent of charge transfer and other 

detailed features of the energy bands. We have found that S is 

between -170 and -70. These results are in general agreement 

with the experimental data on graphite and GICA (nearly universal 

values of £ ~180 relative to TMS and £ v-140). We note that 
x a 

it is very difficult to determine S z of pure graphite experimentally 
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because of the weak signals j tht available experimental data®" 10 ’ 1 ^’ 1 ^ 

are not consistent with each other and further work is needed. 

We will now consider the "dipolar" contributions to CS which 

24 

may be important for GICD. Knight shift measurements indicate 

that there is considerable mixing between the metal conduction 
electrons and the carbon 2px-electrons. The CS would come from 
the dipolar magnetic fields at the nuclear site due to the spin- 
polarized 7T-orbitals in the presence of the external magnetic 
field B. For B parallel to the z-axis, we have 
S z =(^/N) <V*^> 2 p < 3cos 2 0-l> where N is the Avagadro's number, 
yj is the fraction of the polarized electrons that reside on carbon 
atoms, ^ is the Pauli spin suceptibility in cnrVmole and 9 Is the 
azimuthal angle for the <r-orbital of the carbon atom. The first 
pair of angular brackets is radial averaging over the carbon atom 
7T-orbital whereas the second pair is angular averaging. Since 
the angular dependence of the tt -orbital wave function is given 
by cosG, we get 

^3cos 2 0-1> = J c cos 2 9(3cos 2 9-l) sinG d9/ jf eos 2 G sinO dG = 4/5. (14) 

Hence £ z =(4/5) /N) (r'^/^p' Because of the axial symmetry of 
7r-orbitals about the z-axis, it follows that £ =-£/ 2. Thus the 

A Lt 

"dipolar" contribution to CS anisotropy is 

S a = V S x ’ < 6 /5 ) < 7 < /N ) <r _3 > 2p (15) 

where (r~^>2 = 12*2 A" . Because of the close proximity of 

7T-orbital to the carbon nucleus, £ is enhanced by approximately 
a factor of v <r"^> 2p 'v 100 (where v«v 9 t? is the crystal volume per 
carbon atom) over the Pauli contribution to the dimensionless volume 
susceptibility. The experimental susceptibility is the sum of 
several contributions 2 ^, thus the identification of a relatively 
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small Pauli component may be somewhat uncertain. 

We note that the presence of the metallic electrons are also 

1 1 

clearly demonstrated by the electronic specific heat ' which is 
proportional to temperature with the proportionality constant y . 
This constant is nearly proportional to the metal concentration and 
may be related to the Pauli spin pusceptibility ^ i 

^ *(3A*) (pA B ) 2 ir UO 

where |i is the Bohr magnetoni k 0 is the Boltzmann constant, y is 

in units of erg(mole) _1 (°K)’ 2 and ^ is in units of cm^mole" 1 . 

Thus the Pauli susceptibility ^ nay be eliminated between eq. (15) 

and ( 1 6) . The result is £ a *333T^ when S & is expressed in ppm 

and 1 in units of (mJ) (mole) (°K) “ 2 . By combining with the 

experimental value of the "band" contribution S a =-l40 which is 

independent of intercalant concentration, we get 

S a = -140+333^ (17) 

1 1 26 

Using the experimental values ' of y, we have calculated 
S from (17) and have compared with some of the experimental 

cL 

8-16 

values " in the following table for the first and second stage 
GICD with the general composition C X K and C x Cs. From the Knight 
shift data on Cs, it is plausible that the polarized electrons 
are shared equally between the carbon and metal atoms, hence we 
have identified with the mole fraction of carbon atoms in GIC, 
yj =x/(x+l) . The agreement is satisfactory. The variation in 
(isotropic part of CS with TMS as reference) is about a factor of 
10 smaller than the anisotropic part S' • This isotropic term is 
usually ascribed to polarizations of the carbon core s-electrons 
by the unpaired 7T-electron spins , and the factor of 10 reduction 
has been previously observed for similar atoms. 2 ® The results for 
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GICA are also included for comparison. Very small values of y 
have been observed for these compounds, hence comes from the 
"band" contribution only. 



GICA 

CqK 

C 24 k 

CgCs 

C 36 Ca 

V ,mJ mol“*K“ 

2 0 

0.697 

0.241 

0.63 

0.19 

1 


8/9 

24/25 

3/9 

36/37 

9 a , calc. 

-140 

66 

-63 

46 

-79 

$ a , exptl . 

-140 

57 

-63 

35 

-93 

S i( exptl. 

133 

113 

121 

107 

123 


Alternatively, the choice 17= 1 may be inferred from tne x-ray 
29 10 

photoelectron data 7%J which have indicated almost complete charge 
transfer of metallic electrons to carbon atoms. From (17) » the 
calculated & will be changed only slightly (about 20 ppm 

cl 

for CgK and 

agreement with experimental data is still satisfactory. 

In conclusion, the CS in GIC may be expressed as the sum of 
•'band" and’'dipole" contributions. The "band" term is determined 
by the gross features of the two-dimensional band structure and 
is independent of the intercalant concentration. The "dipole" 
term is due to the mixing between metal and carbon electrons and 
is absent for GICA. It is related to the electronic specific 
heat for GICD. There is general agreement with the experimental 
data. 

The pai 4 lal financial support by National Aeronautical and 
Space Administration Grant NAG-5-156 is gratefully acknowledged. 


CgCs, 3 ppm for higher stage compounds) and he 
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